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Abstract. We calculate the neutrino signal resulting from annihilation of secluded dark 
matter in the Sun. In this class of models, dark matter annihilates first into metastable 
mediators, which subsequently decay into Standard Model particles. If the mediators are 
long lived, they will propagate out from the dense solar core before decaying. High energy 
neutrinos undergo absorption in the Sun. In the standard scenario in which neutrinos are 
produced directly in the centre of the Sun, absorption is relevant for E > 100 GeV, resulting 
in a significant suppression of the neutrino spectrum beyond ~ 1 TeV. In the secluded 
dark matter scenario, the neutrino signal is greatly enhanced because neutrinos are injected 
away from the core, at lower density. Since the solar density falls exponentially with radius, 
metastable mediators have a significant effect on the neutrino flux, even for decay lengths 
which are small compared to the solar radius. Moreover, since neutrino detection cross sec- 
tions grow with energy, this enhancement of the high energy region of the neutrino spectrum 
would have a large effect on overall event rates. 
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1 Introduction 

The gravitational capture and subsequent annihilation of dark matter (DM) in the Sun 
offers a compelling possibility for DM detection. DM particles passing through the Sun are 
expected to scatter off nuclei, lose energy, and become trapped by the gravitational field. 
Multiple scatterings cause DM to sink towards the centre of the Sun. If DM self-annihilates, 
the capture is quickly balanced by the annihilation of DM particles. The intensity of the 
annihilation signal is then a probe of the DM scattering cross section on nucleons [1, 2]. 
Signals from DM capture in the Sun (or the Earth) are expected to be detectable in the 
future. Standard Model (SM) particles other than neutrinos, produced in the annihilations, 
interact strongly with the interior of the Sun and are largely absorbed. In this process though, 
they produce high-energy neutrinos which escape, and can be potentially seen by neutrino 
detectors, such as IceCube and DeepCore. Extensive studies of the expected neutrino signal 
from DM annihilation via various SM channels have been performed, in a model-independent 
way [3, 4], and for particular models [5-20]. Limits from Super-Kamiokande and IceCube on 
annihilation of WIMP and Kaluza-Klein DM in the Sun and the Earth have been reported 
in Refs. [21-23]. 

It is possible, however, that DM does not annihilate directly into SM particles, but rather 
into metastable mediators which subsequently decay into SM states, XX ~^ ~^ SM, as 
recently discussed in [24-29]. In such models, the thermal relic WIMP DM scenario can 
be realised as usual, while there is also the potential to explain astrophysical observations, 
e.g. the positron excess observed by PAMELA [26, 27]. The seclusion of DM from the SM, 
which in this class of models communicate only via metastable mediators, can dramatically 
change the annihilation signature of DM captured in the Sun. For example, if the mediators 
are sufficiently long-lived to escape the Sun before decaying, they can produce detectable 
charged-particle or 7-ray fluxes, as discussed in Refs. [30, 31]. 
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If the mediators are short-lived and decay in the interior of the Sun, energetic neutrinos 
remain the only signature. However, even for decay of V inside the Sun, the neutrino signal 
can be dramatically enhanced compared to the standard scenario. This is because high 
energy neutrinos can interact with nuclei and be absorbed before escaping the Sun. In the 
standard scenario, in which neutrinos are produced at the centre of the Sun, the neutrino 
energy spectrum is damped as e"^"/^, with a critical energy of f ~ 100 GeV. However, the 
solar density decreases exponentially with radius, so if neutrinos are injected by V decay 
at larger radii they are subject to much less absorption as they propagate out of the Sun. 
We thus expect the critical energy scale for absorption to increase exponentially with the 
neutrino injection radius, increasing rapidly once the injection point is moved outside the 
dense core. 

For simplicity, we will consider the case where the mediators, V, decay directly to 
neutrinos V — )■ vv', as this will suffice to illustrate the difference between DM annihilations 
with and without mediators. It is of course plausible that the mediators decay to other 
SM model particles, which subsequently produce energetic neutrinos. Although we shall not 
study this case in detail, we note that greatly enhanced signals are also expected in this 
scenario. For instance, if DM annihilates directly to light quarks or muons at the centre 
of the Sun, the neutrino yield is negligible, as these particles quickly lose energy and are 
absorbed before they can decay to produce neutrinos. However, if a mediator decay were to 
inject light quarks or muons at radii beyond the dense core, then a non-negligible neutrino 
flux would result; absorption effects are of course completely avoided if the decay occurs 
outside the solar radius [30]. 

In this paper, we consider the neutrino signal at the Earth from annihilation of secluded 
DM in the Sun, for a range of mediator lifetimes and DM masses. In Sec. 2 we describe the 
production of neutrinos from decay of the mediator, review the relevant aspects of neutrino 
interactions in the Sun, namely absorption, regeneration, and flavour evolution. We also 
describe the numerical approach we follow to calculate the neutrino signal at the Earth. Our 
main results are presented in Figs 4-8. In Sec. 3 we discuss the results, and compare them 
with the standard case of DM annihilation directly into SM particles. 

2 Neutrino production and propagation in the Sun 

During the capture and thermalisation of dark matter in the Sun, the DM particles un- 
dergo multiple scatterings and become concentrated within a rather small region of size 
~ O.Oli?0-y/lOO GeV/m^ around the centre of the Sun, where Rq is the solar radius. The 
finite size of this region has a negligible effect on the final neutrino spectra [1, 3, 4], and thus 
the DM annihilations can be assumed to take place in the centre of the Sun. 
In the scenario considered here, DM annihilates into metastable mediators 



The mediators are emitted with Lorentz factor 7 = m^/my. In the rest frame of the Sun, 
the neutrinos produced from mediator decay uniformly span the energy interval 





which subsequently decay into neutrino pairs 



(2.2) 




(2.3) 
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where (3 is the mediator velocity, and have an angular dispersion 59 ~ I/7 around the 
momentum axis of the mediator [30]. As long as 7 ^ 1, the neutrinos are effectively emitted 
radially outwards, with energies Ey ^ rriy-. This allows us to use 1-dimensional evolution 
equations, similarly to the approach commonly adopted in the standard scenario of DM 
decay directly into SM particles [3, 4, 6]. 

We will focus on DM masses 100 GeV < rriv ^ lOTeV. Neutrinos of energies lower than 
100 GeV have negligible interactions in the interior of the Sun, and their final spectrum is not 
very different from their spectrum at production. Ai Ey>l TeV, the absorption of neutrinos 
in the Sun becomes quite severe, strongly suppressing the signal if neutrinos are injected in 
the centre of the Sun, as is the case for the standard scenario of DM annihilation directly 
into SM particles. However, if DM annihilates into metastable mediators, which travel some 
distance before decaying, the neutrinos transverse smaller optical depth in the Sun and are 
thus absorbed less. As we shall see, this can change the total flux and spectral shape of the 
signal for any DM mass > 100 GeV, and more dramatically so for > 1 TeV. 

Once neutrinos are produced, they undergo charged and neutral current scattering with 
nuclear matter in the Sun, which results in both absorption of high energy neutrinos, and 
re-injection of neutrinos at lower energy. They also undergo flavour oscillations. For the 
energy range of interest, the neutrino flavour oscillations and their scattering interactions 
decouple. The matter potential in the Sun ensures negligible mixing between i/g and v^^Vr 
until the point of the MSW resonance. For 5 GeV Ey < TeV, which encompasses our 
energy range of interest, the resonance occurs after interactions have become unimportant, as 
will be shown in Sec. 2.3. (Neutrinos below 5 GeV make an insignificant contribution to the 
signals we consider.) On the other hand, z^^ and Vr mix maximally, with oscillation length 
typically smaller or comparable to the interaction length. It is then a good approximation to 
assume that oscillations rapidly equidistribute the and Ur fluxes, and thus identify their 
densities and average over their interactions. 

The above considerations allow us to follow the evolution of the flavour eigenstates 
individually, rather than using the full multi-flavour density-matrix formalism. The evolution 
equations have the form 

dr dr 

where p{r, E)dE is the neutrino flux, and j = eor fi, r, for the electron and the averaged muon 
and tau neutrino fluxes respectively. Similar equations hold for the antineutrinos. The flrst 
term on the right-hand side corresponds to the neutrino injection from the mediator decay. 
The second and the third terms describe the neutrino scattering due to neutral-current (NC) 
and charged-current (CC) interactions. Neutrino-neutrino scatterings are unimportant, since 
the neutrino densities are small, and thus the evolution equations remain linear. Each term 
in Eq.(2.4) is described in detail in the following subsections. 

In calculating the neutrino propagation inside the Sun, we adopt the solar nucleon and 
electron density proflles from the Standard Solar Model [32]. For r < 0.9Rq, the nucleon 
density can be well approximated by 

iV5(r) ~ iVoe"^, (2.5) 
with Nq = 1.3 X 10^^ cm~^ and k = 0.1. It will be useful to define the dimensionless optical 



+ 



dr 



+ 



dp^ 
dr 



(2.4) 
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depth (independently of the approximation of Eq. (2.5)): 



1 



o Jo 



Ns{r')dr', 



(2.6) 



which can be inverted to give r{x). For the region of vaUdity of Eq. (2.5), we obtain the 
analytic expressions 



Xa{r) = 1 - e 



2.1 Injection from mediator decay 



ra{x) = kRq In 



1 



1 - X 



(2.7) 



Neutrinos are injected by the decay of the metastable mediators produced in the DM annihi- 
lations in the center of the Sun, as described in Eqs. (2.1), (2.2). In the steady state regime, 
the annihilation balances the DM capture in the Sun, Tann = Cq/2, where the capture rate 
is approximately [2, 5] 



Cq ~ 10^^ s~^ 



100 GeV\ 



(2. 



We will assume that the mediators decay with equal branching ratio to each of the three 
neutrino flavours BRi, = BRg^^^,-. In the rest frame of the Sun, the neutrino energy spectrum 
per V decay is flat, f{E)dE = {BRy/l3my^)dE, for (1 - /3)/2 ^ E/m^ ^ (1 + /3)/2. The 
neutrino injection rate is then 



dp{r, E) 



dr 



CqBR, 1 



/3my /37T 



exp 



Pit 



(2.9) 



We take the lifetime, r, of the mediators to be a free parameter, subject only to the big 
bang nucleosynthesis (BEN) constraint, r < Is [29]. In order to disentangle our results from 
any astrophysical uncertainties and model-dependence, we shall present our results for the 
neutrino fluxes normalized to the product Cq BR,y. 

2.2 Interactions 

Neutrinos in the Sun interact with the nuclei via neutral and charged currents. The interac- 
tion cross sections are rather insensitive to the proton to neutron ratio inside the sun, which 
varies from 2 in the center, to 7 in the outer regions. We thus adopt the constant value 
p/n = 3 throughout. Scattering on electrons is suppressed in comparison to scattering on 
nuclei by the ratio of the electron to nucleon mass, me/mjv, and we will thus ignore it in this 
analysis. 

2.2.1 Neutral-current interactions 

The NC scatterings uiN — t- uiN and i>iN — t- viN shift the neutrino distributions towards 
lower energies. The process can be described as removal of a neutrino from the flux and 
re-injection at a lower energy, as follows 



dp 
dr 



Ns{r) 



/•OD 

-a^^{E)p{r,E)+ dE 
JE 



dE 



p{r,E') 



(2.10) 
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where Ns{r) is the nucleon density of the Sun. The differential and total NC cross sections 
at the energy range of interest have the form [3, 6, 33] 



[a + , (2.11) 

vr 



dE 

a^c{E) = (a + 6/3)^. (2.12) 



NC scattering is flavour-blind, with 

ay ~ 0.06; ~ 0.02; 

(2.13) 

ap = hu, bp = ttu- 
2.2.2 Charged-current interactions 

CC interactions convert a neutrino to an almost collinear charged lepton, viN — )• l~N' and 
viN — )■ l^N'. For electron-type neutrinos, this corresponds to removal of a neutrino from 
the flux. Muon-type neutrinos produce fi^ which decay and re-inject neutrinos in the flux, 
but not before they have thermalised in the plasma and lost most of their energy. The 
re-injected neutrinos with energies ~ few x 10 MeV do not contribute to the high-energy 
neutrino signature of DM annihilations, and for this purpose muon neutrinos that interact 
via CC are considered absorbed. 

Charged-current interactions of Vr and v'j- produce leptons which decay promptly, 
before losing their energy, and re-inject energetic neutrinos in the flux. Besides re-injecting 
a tau- flavour neutrino at a lower energy, the leptonic decays r~ — )• Xi^r (r+ — )• Xv'r) also 
produce a Pe (t'e) 17.8% of the time, and a 1^^ {y^) 17.4% of the time. (The decay 
hadronically with probability 64.8%.) Since the contribution to the i>e i^e) and P^j (v^i) fluxes 
from I'j- (v't-) regeneration is suppressed by the branching fraction ~ 0.18, we will ignore 
this effect^. This stands in agreement with similar neutrino signal calculations in various 
environments which have shown that such a contribution is not significant [34]. 

The CC absorption and regeneration is described by 



dpi 
dr 



= Ns{r) 

CC 



-pi{r,E)a^^^{E)+ / dE' piir, E')fi^i{E, E') 
Je 



(2.14) 



The CC inelastic scattering cross section does not depend on the neutrino flavour, as long 
as we ignore the effect of the r mass. This is a good approximation since the r mass is 
significant only for E < 100 GeV, and at those energies the neutrino interactions in the Sun 
are unimportant. The CC cross section is [3, 6, 33, 35] 

a,,iE) = '-^cE, (2.15) 
vr 

with 

Cy ~ 0.19, Cp ~ 0.13. (2.16) 

For the electron and muon flavours we set fi^i — )■ 0, as discussed above. For tau neutrinos 
these functions encode the appropriate convolution of spectrum arising from CC i'tt'^t 



^Aa with absorption, regeneration has a smaller effect in the mediator scenario compared to the standard 
scenario, so this is an even better approximation tlian usual. 
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interactions, with the regenerated i't^t^t spectrum from the decay. They depend only 
mildly on the incident neutrino energy, and are merely functions of the ratio of the outgoing 
to incident neutrino energy, E/E' . We thus adopt the form of fr-^riE/E') and ff^f{E/E') 
as presented e.g. in Ref. [3] for incident energy E' = 400 GeV. 

As will become clear in Sec. 2.4, the neutrino total interaction strength inside the Sun 
can be conveniently gauged in comparison to the energy scale 



£ = 



'^^>lE^a + b/3 + c)NoKRQ 
vr 



(2.17) 



Using Eqs. (2.13) and (2.16), for neutrinos and antineutrinos respectively, we have 

£^ ~ 140 GeV; £p ~ 213 GeV. (2.18) 

The values of £ in Eq.(2.18) indicate the approximate energy scale at which interactions 
become a significant effect and illustrate that absorption is relevant at lower energies for 
neutrinos than antineutrinos. In terms of £, the neutrino total interaction cross section 
becomes 

The neutrino mean free path. A/, defined by 

r^^' a,ot{E)Ns{r')dr' = 1, 

Jr 

corresponds to optical depth 

6xf=£/E, (2.20) 
and the neutrino interaction probability is 



Pint = ^ ° atot{E)Ns{r')dr' = j{xq - x), (2.21) 



cq ~ 0.9975 is the optical depth at r = Rq. In Fig. 1, we sketch the interaction 
probability contours Pint = 1 and Pint = 0.05, on the E vs r plane. 

2.3 Flavour oscillations 
2.3.1 i^e — i^^i oscillations 

In the inner region of the Sun, the neutrinos experience a matter potential (or refractive 
index) which aligns the flavour eigenstates with (effective) mass eigenstates, and suppresses 
oscillations between electron and mu/tau neutrinos. At lower densities, as the matter poten- 
tial decreases, the e-mu/tau mixing angle increases and the neutrinos go through an MSW 
resonance. The point of resonance occurs when 

A 2 

-^c.V2GFN,{r), (2.22) 

where Ne{r) is the electron density profile of the Sun. The radius at which the resonance 
takes place is shown in Fig. 1. While low energy solar neutrinos go through this resonance 
adiabatically, remaining in a given mass eigenstate, non-adiabatic effects (in which neutrinos 
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Figure 1. Interaction probability contours, Pint — 1 (orange) and Pint = 0.05 (brown), for neutrinos 
(solid lines) and antineutrinos (dashed lines). The black line marks the point of the i'e~'^n,T resonance. 
For the energy range of interest, interactions have effectively ceased when neutrinos undergo resonance. 
We neglect interactions beyond the benchmark point Tc = 0.75Rq (dot-dashed line). Oscillations 
between muon and tau neutrinos occur faster than interactions for energies below the green lines. 
Neutrinos are then equidistributed between the two flavour eigenstates. 




may flip to the other mass eigenstate) are important for neutrinos with energies greater than 
several ~ lOGeV. The level-crossing probability at the point of resonance r = rrcs(-E') is (see 
e.g. [36]) 



PciE) 



o7 cos-^ 6*12 



1 



where 



1 ' 
dlnN, 



dr 



l{E) 



r = rres(-B) 



E ' 



while for antineutrinos 



Pc = Pc [0i2 



12 



(2.23) 
(2.24) 
(2.25) 



The crossing probabilities Pc{E) and Pc{E) are shown in Fig. 2; they become approximately 
constant for E > 100 GeV. 

After resonance, the fluxes of the three mass eigenstates are 



Pi = P»Pc + p!:,1(1-^'c), 
P3 = pSt, 



(2.26) 



where pF(^), p^AE) are the fg and averaged v^-Vt fluxes evaluated at the point of res- 
onance, and similarly for the antineutrinos. We convert the mass-eigenstate into flavour- 
eigenstates fluxes using the standard neutrino mixing matrix [37]. We adopt the best fit 
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Figure 2. Level-crossing probability between mass eigenstates 1 & 2, for neutrinos (solid) and 
antineutrinos (dashed). 



neutrino mixing parameters reported in [37-39], 



Am^i = 7.65 X 10"^ eV^ 



Am 



32 



2.4 X 10"^ eV^ 



sin^ 6^12 = 0.304, 
sin^ 6*23 = 0.5, 
^13 = 0, 
6 = 0. 



(2.27) 



2.3.2 



oscillations 



Muon and tau neutrinos mix maximally, and oscillate with frequency /osc = Am32/2£^. The 
oscillations equidistribute the neutrinos in the two flavour eigenstates, as long as /osc is 
comparable to, or exceeds the interaction rate in the Sun, Tint = crtot{E)Ns{r). In Fig. 1 we 
sketch the maximum energy for which 



/osc ^ Ti 



mt ) 



as a function of the solar radius. This obtains at optical depth 



> 



1 



(2.28) 



(2.29) 



where we used the analytical approximation of Eq. (2.5). Although this condition is not 
strictly satisfied at all depths for the entire energy range of interest, it is easy to show that 
even the highest-energy neutrinos considered here undergo only a small number of scatterings 
< 5, before they lose enough energy and move out in radius to lower solar density, such 
that oscillations become more frequent than their interactions^. Since the mediator decays 



■^On average, every NC interaction in the Sun costs neutrinos about 1/2 of their energy (CC interactions of 
tau neutrinos degrade their energy even more) . After n interactions, the neutrino energy is diminished to En ~ 
Eo/2", and neutrinos have traversed total optical depth Sxtot,n ~ Sxo{2"'^^ — 1), where Sxq = £ /Eq, according 
to Eq. (2.20). If neutrinos are injected in the centre of the Sun, the number of scatterings they undergo till 
they satisfy Eq. (2.29) is found by setting fetot.n = a;oq(£„). This yields n 4.6 + 3.31og(j5o/10 TeV). 
Neutrinos injected away from the centre undergo even less scatterings. 



-8- 



have been chosen to inject equal numbers of all flavours, we don't rely on the oscillations to 
establish the v^-Vr average as an initial condition, but only to maintain it, and thus averaging 
over the muon and tau neutrino populations remains an adequate approximation. 

2.4 Propagation 

We compare the effect of neutrino interactions and oscillations inside the Sun in Fig. 1. For 
all energies in the range 5 GeV ^ < 10 TeV, resonance between the electron and muon 
neutrinos occurs at r > 0.75Rq. At this region, less than 5% of the neutrinos will interact 
before they exit the Sun, Pint < 0.05. This allows for separate treatment of the neutrino 
interactions in the interior of the Sun and the neutrino resonant conversions. 

We calculate the neutrino signal at the Earth according to the following procedure: 

i. We follow the evolution of the neutrino fluxes from the centre of the Sun to Tc = 0.75Rq 
using Eqs. (2.4), taking into account the neutrino interactions as described in Sec. 2.2. 
Eqs. (2.4) are solved analytically as described below. For muon and tau neutrinos 
we set p^^r = = Pj. (and similarly for the antineutrinos), and average over their 
interactions. We neglect interactions at r > Tc, and simply add the contribution to the 
neutrino flux from mediator decays between Tc and the resonance point. 

ii. We evaluate the level crossings at resonance, using Eqs. (2.26). We convert the mass- 
eigenstate fluxes back into flavour-eigenstate ones using the standard mixing parameters 
given in Eq. (2.27). 

iii. We add the neutrino flux generated by decays between the point of resonance and the 
Earth. 

We shall now describe these steps in more detail. 
The evolution equations (2.4) take the form 

/•oo 

-Epj {x,E)+ dE'gj [E/E') p, {x,E') 

JE 



(2.30) 



where £ is given by Eqs. (2.17) and (2.18). For r < Tc, the exponential approximation of 
Eq. (2.5) for the solar profile is valid and the neutrino injection term h{x, E) is expressed 
analytically by 

for the energy interval of Eq. (2.3). The regeneration functions appearing in Eqs. (2.30) are 

ay + hyV? 



9e{u) 

ae{u) 



ttu + K/3 + Cy ' 

ay + hy/?, + Cy 

ay + hyV? 



(2.32) 



ay + 6p/3 + Cp' 

ay + hyV? + Cy/2 ff^fju) 

ay + by/3 + Cy 
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Figure 3. The regeneration functions of Eqs. (2.32), for v^. (purple, solid), (purple, dashed), v^,r 
(blue, solid), P^i.T- (blue, dashed). 

and are sketched in Fig. 3. 

Equation (2.30) can be solved analytically for (7 = and (7 = 1. For a general function 
g, it is possible to solve Eq. (2.30) perturbatively around these two cases. The range of values 
of the regeneration functions, presented in Fig. 3, renders this a suitable treatment for the 
neutrino propagation in the Sun. For the fg and 9e fluxes, g^, ge < 1, and the perturbative 
solution around g = is warranted. For the z/^^r and u^^r fluxes, the perturbative solution 
around g = 1 appears formally more appropriate. We have implemented both solutions and 
checked that, for the DM masses and mediator lifetimes considered here, the differences in 
the final spectra are small and do not qualitatively change our conclusions. For simplicity, 
we will thus adopt the perturbative solution around g = for all fluxes. 

The exact solution of Eq. (2.30) for g = is 

po{x,E)=[ dx' e~^^^T^h{x',E). (2.33) 
Jo 

Let 

p{x,E) = po{x,E) + pi{x,E), 
where pi{x,E) obeys, to 1st order in g, the differential equation 

dpi E 
— = hi{x,E) - —pi{x,E), 

with 

hi{x,E) = -J^ dE' g{E/E')po{x,E'). 
Eq. (2.35) has the same form as Eq. (2.30) for g = 0, and its solution is 

pi{x,E)=[ dx' e-^^^hi{x',E) (2.37) 

Eqs. (2.33), (2.34), (2.36) and (2.37) yield the perturbative solution to the evolution equation, 
to first order in g. 



(2.34) 



(2.35) 
(2.36) 



The contribution to the neutrino flux from decays of the mediators between Tc and 
rresC-B) is 



5p,{E) 



/3my 



exp 



/57T 



exp 



(2.38) 



For 7T > 2 s, the neutrino production beyond the point of resonance is significant and is 
given by 



Cq BRi, 



exp 



exp 



Rse \ 
'hrj 



(2.39) 



with i?sE = 1 A.U. 



The final neutrino fluxes at Earth, calculated as outlined above, are presented in Figs. 4 
- 8 for various choices of the DM mass and mediator lifetime, and discussed in the following 
section. 



3 Neutrino signals 

3.1 Neutrino flux results 

In Figs. 4 - 8 we present the final differential neutrino fluxes at Earth, d^/dE = p{lA.\] ., E), 
calculated as outlined in previous section, for various choices of the DM mass and mediator 
lifetime. All fluxes are normalised to the DM capture rate in the Sun, Cq, and the branching 
ratio BRjy for decay to neutrinos. 

In Fig. 4, we compare the fluxes obtained for different mediator lifetimes, 'jt = 
0.001s, 0.1s, 0.3s, Is, 10s (always assuming 7 ^ 1). For small mediator lifetime, absorption 
effects result in an exponential suppression of the flux at high energy, while down scattering 
or regeneration effects cause a pile up of neutrinos at low energy. These effects can both 
clearly be seen in the left-hand panels of Fig. 4. Relativistic particles transverse the entire 
radius of the Sun in ~ 2.2 s, while the core of the Sun corresponds to ~ 1/10 of this distance 
(cf. Eq. (2.5)). Thus, the smallest mediator lifetime we consider, 7r = 0.001 s, effectively 
corresponds to the standard scenario of DM annihilation directly into SM particles in the 
centre of the Sun'^. For the largest lifetime we consider, 7r = 10 s, almost all the mediators 
decay outside the Sun, so absorption is negligible and the characteristic flat injection spectrum 
of boosted 2-body decays (discussed in Sec. 2.1) is reproduced. Between these two limiting 
cases, the fluxes depend sensitively on the mediator lifetime. 

We can quantify the degree of attenuation in terms of the neutrino injection radius as 
follows: if we inject all neutrinos at radius r = rinj and consider only absorption effects, the 
solution of the flux evolution equation (Eq. (2.30)) is simply 

P = Anje-^^", (3.1) 

where pinj is the flux at the injection radius, and Ax is the optical depth traversed. Neutrinos 
injected at r = rjnj traverse an optical depth of Ax = exp (^—-^;^^ as they propagate to 

^We note that the neutrino spectrum resulting from the prompt decay {'yr < 0.2 s, with 7 S> 1) of each 
metastable mediator to a neutrino pair resembles the neutrino spectrum obtained from DM annihilation to 
relativistic SM particles such as r^r", rather than that for DM annihilation directly to a neutrino pair. In the 
latter case, neutrinos are emitted monoenergetically, which results in a residual peak in the final spectrum, 
at the maximum energy. In the case of decays of boosted particles, neutrinos are injected with an extended 
spectrum, as described in Eq. (2.3). 
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Figure 4. Vf. fluxes for 7T = 0.001 s (solid), 0.1 s (dashed) 0.3 s (dotted) 1 s (dot-dashed), 10 s 
(dot-dot-dashed), and DM masses = 200 GeV (1st row), 500 GeV (2nd row), 1 TeV (3rd row), 
10 TeV (4th row). 
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Figure 5. Neutrino fluxes at tlic Eartli for DM mass = 200 GcV, and mediator lifetimes 7T = 
0.001 s, 0.1 s, 1 s, 10 s. In each graph the lines correspond to: t^e flux (blue, solid), v^^r flux (purple, 
solid), V(, flux (blue, dashed), and P^.,- flux (purple, dashed). 
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Figure 7. Same as in Fig. 5, for DM mass = 1 TeV. 
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r = Rq (see Eq. (2.7)). For = 0, Ax ~ 1, and thus p = p\-a^e £ . More generally, if we 
define a critical energy £' as 



^'(n,j)=£:exp(^-|-J, (3.2) 
the solution for the attenuated flux at the solar surface can be expressed as 

/0 = Anje"£^. (3.3) 

The critical energy £' controls the exponential suppression of the flux at high energies. Since 
£' grows exponentially with injection radius, we see that this critical energy for absorption 
grows rapidly as soon as the injection radius is outside the core, r\a] > kRq. 

In the scenario considered here, mediator decays inject neutrinos not at a single point 
but with a distribution of radii. To account for the absorption we must employ the injection 
function of Eq. (2.9), or equivalently 



exp -- — exp --pq-r-; dr. (3.4) 



This provides a good approximation for the flux for high energies E > £, where absorption 
is the leading effect. Regeneration effects become relevant at lower energies, where a pile up 
of low energy neutrinos enhances the fluxes above the levels indicated by Eq. (3.4). 

Annihilation via the mediators enhances the prospects for detection of a DM annihilation 
signal from the Sun. In the right-hand panels of Fig. 4 (and in Figs. 5 - 8) we plot the quantity 
Ed^/dE. As the neutrino scattering cross-section on nucleons is proportional to the energy 
(see Eq. (2.19)), the area under the Ed^/dE vs. E graphs is a guide to the relative event 
rates in neutrino detector. For all DM masses and mediator lifetimes, the prospects for 
detection are better at neutrino energies E^, > 100 GeV, and dramatically so as the DM 
mass is increased. Since the high energy neutrinos are the most enhanced and also have the 
largest detection cross sections, the mediator scenario is particularly effective in enhancing 
event rates. 

The neutrino spectral shape is a striking signature of DM annihilation via metastable 
states. In the standard scenario, there is never significant neutrino fiux at energies Ei, > 
1 TeV. If DM annihilates directly into SM particles, Ed^/dE exhibits a peak around ~ 
100 GeV (see Figs. 5-8). The peak occurs because the absorption increases strongly with 
energy and always quells the spectrum at higher energies. In fact, for > 1 TeV, neutrino 
interactions in the Sun drive the spectrum to a fixed shape independent of the DM mass 
or annihilation channel (dubbed the "limit spectrum" in Ref. [3]) with a sharp cutoff at 
~ 1 TeV. These features are reproduced by our calculations for 7T = 0.001 s. In the 
presence of mediators, however, the peak of the spectrum moves to higher energy if 7T is 
small, and completely disappears for 7T > 0.3s (which corresponds to mediator decay outside 
the core of the Sun). Although the small shift in the spectral shape for small 7r may not be 
readily discernible, the rising trend of the spectrum [Ed^ / dE) at high energies that occurs 
in the case of sufficiently long-lived mediators is a firm signature of this scenario. It cannot 
arise in the standard scenario, even for very large DM mass. A potential observation of a 
neutrino signal from the Sun at energies Ey>^\ TeV can thus be justified only in the context 
of DM annihilation into metastable states with 7r > 0.1 s. 

We compare the fluxes of all neutrinos flavours, for both neutrinos and antineutrinos, 
in Figs. 5-8. We have assumed that the mediators decay with equal branching ratio to 
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all flavours. However, interactions induce flux differences between neutrinos and antineutri- 
nos, and between the different flavours. Neutrinos interact more strongly in the Sun than 
antineutrinos, thus are absorbed more, resulting in smaller fluxes. Due to the Vr and Vr re- 
generation, the muon/tau neutrino and antineutrino fluxes are higher than the electron- type 
ones. These features are more apparent for larger DM masses and smaller mediator lifetimes, 
when interactions are stronger. As the mediator lifetime is increased, such that interactions 
become less important, the spectra for all flavours of v and v merge. 

3.2 Muon event rates 

We shall now estimate the number of muon track events in the IceCube detector for a range of 
mediator lifetimes, relative to the scenario in which the neutrinos are produced in the solar 
core. We estimate event rates following the procedure in Refs. [40-42], which is outlined 
below. 

A muon neutrino of energy E produces a muon with energy = E[l — y) when it un- 
dergoes a charged current interaction, where y is the charged current inelasticity parameter. 
For simplicity we will adopt the mean values of y tabulated in Ref. [40], which are approx- 
imately ~ 0.45 for neutrinos and ~ 0.35 for antineutrinos. Muons propagate over a long 
distance before they decay, such that IceCube can detect muons created from interactions 
well outside the detector. The muons lose energy as they propagate, and have a range given 
by 



R,{E,,E'^^) = -^\n 



a + pE^ 

thr 
ithr 



a + pE^ 



(3.5) 



where a = 2.0 MeV cm /g, /3 = 4.2 X lO"'' cmVg, and El^' is the muon energy detection 
threshold for the detector (see, e.g., Refs. [40-43]). We wih take E]^"" to be 100 GeV. The 
probability that a neutrino of energy E makes a muon track that is detected with energy 
above £^1','^'^ is then 

PiE, E'l^') = pjviVAa„ucleoni?M(^(l " y)^^^'), (3.6) 

where pat is the target nucleon density and Na is Avogadro's number. The muon track event 
rate is thus proportional to 

^^j^P{E,Ef^')A''^dE (3.7) 
dE 

where A*^^ is the effective area of the detector, which is ~1 km^ for IceCube [44]. 

In table 1, we show the relative number of muon track events for different choices of 
the mediator lifetime, normalised to the case where r ~ (i.e., production in the centre 
of the Sun). For the latter, we adopt the fluxes obtained here for 7r = 0.001 s, in order to 
retain consistency with our calculations. Both neutrino and antineutrino events are included, 
weighted according to their respective CC interaction probabilities (which have a ratio of 
approximately 2:1 in ice). For dark matter masses close the transparency energy of the 
Sun ~ 200 GeV, the effect of the metastable mediator is to enhance the flux by about a 
factor 2. For larger dark matter masses, well above the energy at which the Sun becomes 
opaque neutrinos, annihilation via metastable mediators has a huge effect on the flux, with 
an enhancement of an order of magnitude for lifetime as short as 0.1s, and an enhancement 
of about 1000 for lifetime comparable to the solar radius. Particularly for these heavy dark 
matter masses, the presence of the metastable mediator greatly enhances the flux which 
emerges from the Sun, and therefore the detection prospects. 
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500 GeV 


1 TeV 
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0.3 s 
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7.6 


2.7 X 10^ 


Is 


2.3 


4.8 


1.4 X 10 


1.4 X 10^ 


10 s 


2.9 


6.7 


2.2 X 10 


3.6 X 10^ 



Table 1. The ratio of events observable in the presence of mediators of lifetime 7r over those in 
the absence of metastable mediators, $(7r)/$(7T — >■ 0), integrated over muon energies > lOOGeV, 
for various dark-matter masses and mediator lifetimes. We have used the calculated fluxes from our 
analysis for = 0.001 s as the benchmark fluxes for DM annihilation directly into SM particles. 

4 Conclusions 

The search for high energy neutrinos produced by the annihilation of DM captured in the 
Sun is a promising means of indirect DM detection. A future detection can provide infor- 
mation about the DM annihilation channels, which is critical for the identification of the 
particle nature of DM. It is possible, however, that DM does not annihilate directly into SM 
particles. We have examined a scenario in which DM annihilates to metastable mediators, 
which subsequently decay to SM states. In these models, though the DM and SM sectors are 
secluded from each other and communicate only via the mediators, the standard thermal- 
WIMP scenario can be realised as usual. Existing DM indirect detection bounds are modified 
in these models. Particularly interesting is the signal from DM capture in the Sun, which 
is enhanced at high energy and provides distinct features which can unmask the secluded 
nature of DM. 

In the case where the mediators propagate out from the dense solar core before decaying, 
the neutrinos produced in the decays experience a much smaller optical depth and thus much 
less absorption, compared to neutrino production at the centre of the Sun. In the standard 
scenario, the absorption is so severe that the flux is almost completely suppressed beyond 
about 1 TeV. The observation of a solar neutrino flux beyond 1 TeV would thus be striking 
indication of the presence of metastable mediators. These features would be discernible even 
with small number of events, or poor detector energy resolution, and remain valid independent 
of the specific decay channel of the mediators into SM states. Finally, this enhancement of 
the neutrino flux at high energy is particularly effective in improving detection prospects, 
given that the detection cross sections grow with energy. 
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